Phase separation in mixtures of colloids and long ideal polymer coils 
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Abstract 

Colloidal suspensions with free polymer coils which are larger than the colloidal particles are considered. 
The polymer-colloid interaction is modeled by an extension of the Asakura-Oosawa model. Phase separation 
occurs into dilute and dense fluid phases of colloidal particles when polymer is added. The critical density 
of this transition tends to zero as the size of the polymer coils diverges. 



By the addition of nonadsorbing polymer, colloidal 
suspensions can be made to separate into a dilute 
phase, a colloidal 'vapour', and a dense phase, a col- 
loidal 'liquid' [fflj. Nonadsorbing polymer does not 
adsorb onto the surfaces of the particles, the inter- 
action between a monomer and colloidal particle is 
repulsive. The colloidal particles are spheres of di- 
ameter a, and we can characterise the size of a poly- 
mer coil with its root-mean-square end-to-end separa- 
tion, Re- For polymer coils smaller than the particles, 
Re < c, the polymer-colloid interaction is, to a not 
unreasonable first approximation, a hard-sphere re- 
pulsion: the Asakura-Oosawa model |2], [|. For values 
of Re/c = O(0. 1), it has been found that to a good 
approximation the effect of the polymer is to create 
a short-ranged, range 1Z ~ a + Re, pairwise additive 
attraction between the colloidal particles |§,]£]]. If this 
range is not too short it induces a phase separation 
into dilute and dense fluid phases of the colloidal par- 
ticles. In the other limit, that of large values of the 
ratio Re/ct, it is known that the polymer creates an 
effective attraction between the colloidal particles but 
that this is not pairwise additive |,g. Here, we ex- 
tend the Asakura-Oosawa model to deal with values of 
Re I® > 1) and then go on to show that the phase be- 
haviour for large Re/& is very different from that in- 
duced by a pairwise additive attraction. As Re — > oo 
the colloid density at the critical point tends to zero, 
whereas with a pairwise additive attraction, however 
long the range of the attraction 1Z, the density at the 
critical point remains nonzero. A van der Waals fluid 



has a pairwise additive attraction of infinite range but 
a critical volume fraction of 0.13 [J5J. The colloidal 
particles can be smaller than the polymer coils either 
because the particles are small, only a few nanometers 
across, the typical size of a globular protein ||, or the 
polymer coils are large, as they are for DNA ||. 

The colloidal particles are modeled by hard 
spheres, there are no attractions between them. The 
interaction potential between two colloids, ucc, is 
then 
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where r is the distance between their centres. The 
polymer coils are taken to be ideal and so do not in- 
teract with each other, i.e., upp(r) = 0. For the in- 
teraction between a colloidal particle and a polymer 
coil, ucp, we start from the Asakura-Oosawa model 
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where op is an effective diameter of the polymer coil; 
it is close to Re, the root-mean-square end-to-end sep- 
aration of the polymer. The polymer does not adsorb 
onto the particles; the interaction between a monomer 
and the surface of a particle is repulsive. When the 
polymer is no larger than the colloid the Asakura- 
Oosawa model is reasonable. In the limit a 3> Re, 
the colloidal particle resembles a hard planar wall on 
the relevant length scale for the polymer, Re, and a 
hard wall excludes a polymer from a slab of height of 



order R B [1C, 11]. However, in the opposite limit, that 
of Re 3> c, the Asakura-Oosawa model is incorrect, it 
predicts that a particle excludes a polymer coil from a 
volume of R E , whereas if Re is much larger than the 
diameter of the colloid the colloid-polymer interaction 
must be extensive in the length of the polymer. For 
an ideal polymer Re = an 1 / 2 , where a is the monomer 
length and n is the number of monomers. The inter- 
action must be extensive in n and so scales as R E not 
R E pO|,[Tl[|. Thus, we cannot use the Asakura-Oosawa 
model for long polymers. We propose an extended 
Asakura-Oosawa model to deal with the case Re > cr. 
Figure |l] is a schematic of the proposed model. 

First we rescale the monomer size to the colloid 



diameter a [12]. As both Re and the exponent of one 
half remain constant when we rescale the monomer 



size, we have R\ 
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a rib, where rib is the num- 



ber of blobs: effective monomers of length a. This 
yields 



n b = Re/o 2 , n B > 1. 
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So, we now have a polymer of rib blobs, each of which 
is a across. Each of these blobs is no larger than the 
colloidal particle so the Asakura-Oosawa interaction, 
Eq. (|2j), is a reasonable (although not a quantitative) 
description of the interaction of a single blob with a 
colloid. Thus, our model for a long polymer is an ideal 
chain of blobs of diameter a, each of which interacts 
with a colloidal particle with an interaction potential 
given by Eq. (g) with cr p = a. 

For ideal polymers, when calculating phase dia- 
grams, it is simplest to work in a semigrand ensemble 
fl |J, |bi. This is an ensemble in which the variables 
are the volume, V, the number of colloidal particles, 
Nq, and the activity of the polymer, z. As all our 
interactions are athermal, the temperature T is a not 
a relevant variable. For simplicity we will use units 
such that the thermal energy kT = 1. The number 
density of the colloidal particles pc = Nc/V. We 
use a reduced density for the colloidal particles, the 
volume fraction r/ = (Tr/6)(Nc/V)a 3 , and a reduced 
activity of the polymer z* = zR E . We also define a 
reduced density of polymer, p* = ppR E , where pp is 
the number density of polymer coils. The colloidal 
volume fraction is close to one half when the colloid 
(in the absence of polymer) crystallises, and p* is close 
to one when the polymer coils start to overlap, and is 
less than one in the dilute regime Jj~2[ ]. 

The semigrand potential per colloidal particle uj is 
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where a B s is the Helmholtz free energy per particle 
of hard spheres, and a(w;n B ) = exp [ — fJ-ExiVi n B)], 
with pEx(v> n B) the excess chemical potential of a 
chain of iib blobs in a system of hard spheres at a 
volume fraction r\. Equation (0) for the semigrand po- 
tential is approximate, it is essentially the semigrand 
potential expressed as the z = limit (no polymer) 
plus a series expansion in powers of z, truncated after 
the first term. Thus, it becomes less accurate as the 
polymer activity increases. It can be derived in a cou- 
ple of ways, see Refs. |, g, H U, . The Carnahan- 
Starling equation of state |15| is known to be accurate 
so we will use thermodynamic functions and correla- 
tion functions derived from it. The density of polymer 
coils when its activity is z* is 



p* = z*a(n,n B ), 
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which is just the definition of the excess chemical po- 
tential, rearranged. 

If riB = 1, then the polymer-colloid interaction is 
modeled by just one hard sphere of diameter cr, as in 
the original Asakura-Oosawa model. Then p B x is just 
the excess chemical potential of hard spheres, which 
may be easily derived from the Carnahan-Starling 
equation of state. For n B > 1 we require the excess 
chemical potential of a chain of blobs. Fortunately 
the problem of calculating this quantity has occurred 
in the treatment of dense liquids of oligomers, such as 
alkanes, and polymer melts. One of the best known 
theories is that of Wertheim [|0|, which he termed 
thermodynamic perturbation theory 1 (TPT1). Its 
prediction for pEx{ili n B) is p^-p^1 



(HS) , 
PEX = n BP E X ~ \ n B 



where p^x* anc ^ 9hs, are the excess chemical poten- 
tial and pair distribution function at contact, respec- 
tively, of hard spheres, pex is the work done in insert- 
ing a chain of n B spheres of diameter a into the fluid 
of hard spheres, which is equal to the work done in 
inserting n B widely separated spheres (the first term 
on the right hand side of Eq. (|6|) plus the work done 
in bringing the n B spheres together into a linear chain 
of spheres at contact (the second term) |lS|j . 

Equations (0) and (^) are all that is required to 
calculate the phase diagram in the rj — z* plane. Then 
Eq. (||) can be used to calculate the polymer den- 
sity from its activity and so these diagrams may be 
mapped onto the rj — p* plane. The pressure and the 
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chemical potential of the colloid are both derivatives 
of Q, and so may be determined from Eq. (Q), and 
then the conditions of equal chemical potential and 
pressure may be solved for phase coexistence. Results 
for us = 5 are shown in Fig. |2[ This corresponds to 
Re = y/^cr, the end-to-end separation of the polymer 
is a little over twice the diameter of the colloid. If the 
colloid is a globular protein molecule this corresponds 
to an Re of around lOnm. 

The density of polymer coils decreases exponen- 
tially with increasing colloid volume fraction, with a 
coefficient in the exponential which is linear in rig, Eq. 
(|6|). Thus for polymers which are several blobs long, 
the polymer density at a polymer activity z* = O(l) 
and at high colloid volume fractions, rj > 0.3, is ex- 
tremely small. Thus, crystallisation of the colloid 
takes place in the presence of almost no polymer and 
so occurs very close to its value for hard spheres, which 
is 77 = 0.49 m, so is almost completely unaffected at 
these polymer activities. We do not show the crys- 
tallisation transition in Fig. || because it is at much 
higher densities than the vapour-liquid transition but 
the density of the fluid phase which coexists with the 
crystalline phase is essentially a vertical line at a vol- 
ume fraction of 0.49. 

In Fig. [2] the volume fraction of the colloid at the 
critical point is very low. The density range of the 
colloidal liquid is very large, from the volume frac- 
tion at the critical point, 0.048, to 0.49. In Fig. | 
we have plotted the volume fraction of the colloid 
at the critical point, r]cp, as a function of polymer 
size, tib- For large ng, it decreases as n^ 1 (i?^ 2 ). 
Thus it tends to zero as Re — * 00 unlike the case 
for a pairwise additive attraction where as its range 
1Z — > 00, the volume fraction at the critical point tends 
to 0.13. The mixture phase separates into colloid- 
rich (polymer-poor) and colloid-poor (polymer-rich) 
phases at very low colloid concentrations when the 
polymer is larger than the colloid. The reason for the 
n^ 1 scaling is clear from Eqs. (|5|) and ([]). The phase 
with the higher colloid density, the colloidal liquid, 
must have a sufficiently high density that at fixed z* 
the polymer density is significantly below that in col- 
loidal vapour. Now, from Eq. (||) we see that this 
requires a [lex which is at least of order unity (re- 
call that kT = 1) in the liquid phase. At low colloid 
density /iex = (A + Bub)^ + 0(rj 2 ), where a and b 
are constants. Thus, the colloidal volume fraction at 
which fiEX = 1 varies as I/^A+Btib), which gives rise 
to a critical density with the same scaling. A critical 
volume fraction scaling as R~^ 2 is consistent with work 
on a polymer molecule in the presence of a density of 



fixed obstacles [20| which finds that the reduction in 
entropy of the polymer molecule is of order one when 
the number density of obstacles of diameter a is of 
order l/(aR 2 E ). 

A further point to note is that as the activity of 
the polymer is increased the density of the colloid- 
rich phase increases rather slowly . In Fig. ^, even 
when the polymer activity is twice that at the critical 
point, the colloid volume fraction in the colloid-rich 
phase is only around 0.15. Ultimately, we expect that 
if the polymer activity is high enough there will be a 
triple point, where the dense fluid is sufficiently dense 
that it coexists not only with a dilute fluid phase but 
with a crystalline phase. However, this will be for 
much larger polymer activities than shown in Fig. y. 
The polymer density in the colloid-poor, polymer-rich 
phase will be many times the overlap concentration, 
p* = 1. Although simultaneous coexistence of dilute 
and dense fluid phases and a crystalline phase have 
been observed in experiment for colloid + polymer 
mixtures p], this has been for values of Re no larger 
than the colloid diameter. Observing simultaneous co- 
existence of these three phases in experiment may be 
difficult if Re is significantly greater than the colloid 
diameter. 

In the limit of short polymers, Re a few tenths 
of a or less, the effect of polymer is to induce an at- 
traction which is short ranged, 1Z ~ a + Re, and to a 
good approximation pairwise additive || . If the phase 
diagrams for colloid + short polymer ||,13], and for 
particles with a short-ranged pairwise-additive attrac- 
tion p2| , are compared they are seen to be qualita- 
tively the same. In particular, in both cases as Re 
or 1Z shrinks, fluid-fluid coexistence disappears from 
the equilibrium phase diagram. Thus, assuming that 
free polymer is equivalent in effect to a pair attraction 
between colloidal particles is a reasonable assumption 
for small values of Re/cf but not for large values. 

We have proposed an extended Asakura-Oosawa 
model to model the interaction between colloidal par- 
ticles and ideal nonadsorbing polymer coils with end- 
to-end separations Re larger than the diameter of the 
colloidal particle a. For a globular protein with diam- 
eter of a few nms, this would mean a polymer molecule 
with an Re of 5nm or more. As with smaller poly- 



mer coils [13 1, the polymer induces a vapour-liquid- 
like separation into two fluid phases. One rich in col- 
loid (a colloidal 'liquid') but poor in polymer and one 
poor in colloid (a colloidal 'gas') but rich in polymer. 
We showed that the critical point of this transition 
moves to lower and lower colloid densities as the poly- 
mer coils become larger. This is qualitatively different 



from what would be found if the effect of polymer was 
to induce a pairwise additive attraction between the 
colloidal particles. Let us compare the variation of 
the colloid density at a critical point induced by a 
pairwise additive attraction of varying range 1Z, with 
that induced by a polymer of varying size, Re- We 
find that for small 1Z or Re the variation of the criti- 
cal density is similar in both cases. The critical den- 
|23| as 1Z or Re shrinks; if not 



sity increases [13 



preempted by crystallisation it tends to the random- 
close-packed density of hard spheres as 1Z or Re — ► 
|23| . However, in the other limit, that of large 1Z or 
Re, the variation in the critical density is very differ- 
ent in the two cases. With a pairwise additive poten- 
tial it shrinks to 0.13 [0 as 1Z — > oo and then goes 
no lower, whereas with polymer coils, the present the- 
ory predicts that the critical density tends to zero as 
Re — ► oo. 
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Figure 1: A schematic of the extended Asakura-Oosawa model. The black discs represent the colloids and 
the curve represents a polymer coil. The rescaled monomers used to estimate the polymer-colloid interaction 
are drawn as dashed circles. 




Figure 2: The phase diagram of a colloid + polymer mixture, with a polymer of ub = 5 blobs of diameter 
equal to that of the colloid, a) is the diagram in the rj — z* plane, and b) is it in the rj — p* plane. The dashed 
lines are tie lines connecting coexisting phases and the circle in b) marks the critical point. 
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Figure 3: The volume fraction of colloidal particles at the fluid-fluid critical point, r)cp, as a function of the 
size of the polymer, measured by ng- 
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